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Abstract: We solve the Tolman-Oppenheimer-Volkoff equation using an equation of state 
(EoS) calculated in holographic QCD. The aim is to use compact astrophysical objects like 
neutron stars as an indicator to test holographic equations of state. We first try an EoS 
from a dense D4/D8/D8 model. In this case, however, we could not find a stable compact 
star, a star satisfying pressure-zero condition with a radius R, p(R) = 0, within a reasonable 
value of the radius. This means that the EoS from the D4/D8/D8 model may not support 
any stable compact stars or may support one whose radius is very large. This might be due 
to a deficit of attractive force from a scalar field or two-pion exchange in the D4/D8/D8 
model. Then, we consider D4/D6 type models with different number of quark flavors, 
Nf = 1,2,3. Though the mass and radius of a holographic star is larger than those of 
normal neutron stars, the D4/D6 type EoS renders a stable compact star. 
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1. Introduction 

The approaches based on the Anti de Sitter/conformal field theory (AdS/CFT) correspon- 
dence [1, 2, 3] are showing many interesting possibilities to explore strongly interacting 
systems such as dense baryonic matter, stable/unstable nuclei, and strongly interacting 
quark gluon plasma. The basic strategy is to introduce an extra space, which roughly 
corresponds to the energy scale of 4D boundary field theory, and try to construct a 5D 
holographic dual model that captures certain non-perturbative aspects of strongly cou- 
pled field theory. There are in general two different routes to modeling holographic dual 
of quantum chromodynamics (QCD). One way is a top-down approach based on stringy 
D-brane configurations. The other way is so-called a bottom-up approach, in which a 5D 
holographic dual is constructed from a boundary field theory such as QCD. For reviews on 
holographic QCD we refer to [4]. 

Neutron stars are the most important astrophysical compact objects in which the dense 
matter physics can be explored and tested. Since one cannot obtain the information on 
the inner structure of neutron star directly from observation, mainly masses and radii of 
neutron stars are being considered. Due to the uncertainties in dense matter physics at high 
densities, various equations of states have been suggested [5]. One of the main problems 
lies in the fact that most of the microscopic theories are based on the data at low densities, 
only up to normal nuclear matter density. Hence, in order to guide high density behavior, 
new approaches which are based on QCD fundamental symmetries have been suggested. 
This work is along this line of work even though it is still far from realistic. 

Most accurate measurements on the neutron star masses has been done in double 
neutron star binaries and well measured neutron star masses in these systems are consistent 
with < 1.5Mq [5]. However, recently, a 1.97 ± 0.04M Q neutron star has been observed in 
a neutron star - white dwarf binary J1614— 2230 using Shapiro delay [6]. In addition, both 
masses and radii of neutron stars in low mass X-ray binaries have been estimated and more 
data will come soon [7, 8]. These observations already ruled out some sets of equations 
of state. In the future, once the advanced LIGO is operating, gravitational waves from 
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neutron star binary inspirals are expected to be detected. These observations may be able 
to provide clues on the inner structure of neutron stars. 

In this work, to see what the equation of state (EoS) from holography says about 
nature, we solve the Tolman-Oppenheimer-Volkoff (TOV) equation using the EoS derived 
from a top-down model. We use compact astrophysical objects like neutron stars as a 
testing agent. We first take an EoS from a cold and dense D4/D8/D8 model. With this 
EoS, we cannot find any stable compact stars with a reasonable radius. This might be 
due to a weakness of attractive force in the D4/D8/D8 [9]. In nuclear physics, the long 
range attractive force is mostly given by one-pion exchange, while the intermediate range 
~ 1 fm attraction is dominated by a scalar field or two-pion exchange. Note that in the 
D4/D8/D8 [9] a scalar that couples to baryon fields is missing. Then, we consider a D4/D6 
type model with different number of quark flavors, Nf = 1,2,3. Since the D4/D6 type 
models have both vector and scalar bulk fields, we might expect that the EoS calculated in 
this model could support a stable compact star. We find that the D4/D6 type EoS allows 
a stable compact star, though the mass and radius is a bit larger than observed neutron 
stars. 

One cautionary remark is that a simple classical analysis in the holographic dual model, 
both top-down and bottom- up, captures only the leading N c contributions, and we are 
bound to suffer from sub-leading corrections. For instance for low densities we may neglect 
a back-reaction from a bulk U(l) gauge field whose boundary value defines a chemical 
potential and corresponding density. However, at high densities, though we have no clear- 
cut number to tell low density from high density, the back-reaction will not be negligible. 

2. D4/D6 model 

We briefly summarize the D4/D6 system [10] and comment on its extended version in 
dense matter with Nf > 1 [11, 12]. The model contains N c number of D4 branes and Nf 
flavor D6 branes whose configurations are given in Table 1. In the probe limit, the N c D4 
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Table 1: The brane configurations : the background D4 and the probe D6 



branes are replaced by corresponding type IIA SUGRA background, while Nf D6 branes 
are treated as probes. The solution of confining D4 brane reads 



= Q 3/ \^ dx v + f(u)dT>) + (§) 3/2 (^ + u*dn 
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where f(U) = 1 — (Ukk/U) s . The D4 branes are compactified on the r-direction, and U is 
considered as a radial coordinate of the transverse to the world- volume of branes. 64 and 
f^4 denote the volume form and volume of the unit £ 4 , respectively. For convenience, we 
introduce the dimensionless radial coordinate £ defined as 

The parameters in string theory are related to those of the gauge side as 

L ~ 2M^ ' UKK ~ 9 ' 9s ~ 2nN c M KK l s ' {2A) 

with the 'tHooft coupling A = gy M N c . To describe the D6 brane configuration, the back- 
ground metric (2.1) is rewritten as 



ds 2 



j) 3/ \v^dxV + f(U)dT 2 ) + ^y /2 (^y(dz 2 + z 2 dnl + dr 2 + r 2 d4> 2 ) (2.5) 



where £ 2 = z 2 + r 2 . The location of D6 brane is described by r(z) with = and r = tq. 
The distance between D4 and D6 branes is interpreted as a quark mass; r(oo) ~ m q . Then 
the induced metric on the D6 brane is 

dS ° 6= {jf /2 V^dxV+ ^y /2 ^Ly^l + r' 2 )dz 2 + z 2 dn 2 ] , (2.6) 

where the prime denotes the derivative with respect to z. The embedding configuration of 
probe brane is governed by Dirac-Born-Infeld (DBI) action. The DBI action for D6 brane, 
in the case of Nf = 1, becomes 

S D 6 = -W>y d 7 a e _ V-det(P[g] + 2to'F) 

= -re J dtdz z 2 (l + l/£ 3 ) 4/ y (1 + l/£ 3 ) 4/3 U + r '2) _ F i ( 2 . 7 ) 

where ^ 1 = (2vr)% 7 , r 6 = ^gJ^V^U^/ '4) and F = (2 2 / 3 /U KK )27ra' 'F zt = A' t . Here, 
for later convenience, more general DBI action including a gauge field At is considered. 
The source of the gauge field is the end points of fundamental strings. We define the 
dimensionless quantity Q, 



dF 



1 + l/£ 3 ) 4 / 3 (l + r' ) — F 2 



and Q is related to the number of fundamental strings Q by Q = Q / {2hoi't§)(Ukk/2 2 ^) . 
Here Q is identified with the total quark number of a system. The Hamiltonian of the D6 
brane is obtained through the Legendre transformation. 

H m = r e Jdz ^(1 + l/£ 3 ) 4 / 3 [z 4 (l + 1/e 3 ) 8 / 3 + Q 2 ] (1 + r' 2 ) (2.9) 
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This D4/D6 model is extended to study dense matter in confined phase [11] and 
to investigate a model for transition from nuclear matter to strange matter [12]. Here 
we briefly review the D4/D6 model in dense matter and refer to [11, 12, 13] for details. 
In holographic QCD, a compact D4 brane wrapping on the 4-sphere 5 4 transverse to 
IR 1,3 is introduced as a baryon [13]. To describe the compact D4 brane, the S 4 part of 
the background metric (2.1) is rewritten as dO| = d9 2 + sin 2 0df2| . The world- volume 
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Table 2: The brane configuration : the compact D4 



coordinate for the compact D4 is ( t, 8, (p a ) where is the polar angle from the south pole. 
To study dense matter, we turn on the time component of the gauge field At. The 50(4) 
symmetry gives the ansatz that all fields depend on 9, that is, £(#) and At(6). The induced 
metric on this D4 brane is 

ds 2 m = -(j) 3/2 dt 2 + (L 3 Uy/ 2 [(l + pd9 2 + S m 2 9dn 2 ] (2.10) 

where £ = d^/dO. The DBI action of the compact D4 brane is 

5 D4 = -fi 4 J d 5 a e - V-det(P[s] + 2ira'F) 

= -T 4 J dtdO Sm 3 9^J (1 + 1/^3)4/3(^2 + £2) _ p2 ( 2 .11) 

where ^ = (2vr) 4 / s 5 , r 4 = // 4 57 1 ^3^ 3 (^kk/2 2/3 ) and F = (2 2 / 3 /U KK )27ra' F et = A t . 
The dimensionless displacement is defined as 

9C <» Sm3$F -D, (2.12) 



which is given by Chern-Simons (CS) term. 

S^ 4 S) = A l f\G i = t 4 J dtdO 3 sin 3 # A t (2.13) 

From the equation of motion of the gauge field, D is obtained explicitly. 

D{6) = 2- (3 cos 0- cos 3 0) (2.14) 
The Hamiltonian of the compact D4 brane is also given by the Legendre transformation. 

H m =T A Jd6 + 1/^3)4/3(^2 + ^2)^2 + gin 6 Q 



We consider the system in which the number of fundamental strings is Q. Since N t 
fundamental strings are attached on each compact D4 brane, there exist Nb(= Q/N c 
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number of D4 branes. The other end of the fundamental string is located on the D6 brane 
as a source of the gauge field. Since the tension of the string is always larger than that of 
the brane, the fundamental strings pull both D4 and D6 branes. Eventually the length of 
the string vanishes and the compact D4 and D6 branes meet at a point. To be a stable 
system, the force at the cusp should be balanced: (Q/-W c )/d4 = fm[Q] • Force due to the 
brane tension at the point is given by the variation of its Hamiltonian with respect to the 
position of contact point. 

f=%T ( 2 - 16 ) 

OC, c fixed other values 

At zero temperature, energy density e and pressure p are calculated as 

e=^(§-n m + n r °llQ}) and p=l^ e |[Q] (2.17) 

where and Cl denote the grand canonical potential and its Legendre transformed one, 
respectively. These quantities can be obtained as on-shell action and Hamiltonian. 

& = S'lon-shell , ^ = -^Ion-shell (2-18) 

When we calculate the on-shell values of probe D6 branes, the integration result with the 
infinite range of z diverges. Therefore, we should subtract out the zero-density contribution. 

rcg - [Q] = n[Q\- n[Q = 0] (2.19) 

The baryon number density is also given by p = Q/(N c Vs) . 

If we want to extend the single flavor description to a multi-flavor case (Nf > 1), it is 
convenient to introduce the quantity q = Q/Nf . In this work we consider all Nf D6 probe 
branes have identical asymptotic heights, i.e., Nf quark flavors with the same masses. In 
this case we can trivially extend the single flavor case to the multi-flavor one. For a fixed 
Q, Q/Nf number of fundamental strings are attached on each D6 probe brane. Then the 
force balancing condition is modified as (Q / 7V C ) /d4 = Nff^lq] . The energy density and 
the pressure are 

( — tin* + Nr » \a]\ . v=^- 



i xv n m + Nf^f[ q ]j , p =-L^I[ q }. (2.20) 
3. Numerical Results 

Now, we use the EoS obtained in holographic QCD models to calculate the mass-radius 
relation of a compact star. To this end, we solve the general-relativistic TOV equation 
using the holographic EoS: 

dp 1 2m + 8vrr 3 p 

Tr=-2 {e + P) r(r-2m) ' ^ 

where r is the circumferential radius, and geometric units, c = G = 1 are used. The energy 
density, e, is given by e = (E/A + m^p, where EjA is the energy per baryon, = 938 
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MeV is the mass of the baryon, and p is the baryonic number density. The integrated 
mass of the compact star is given by m(r) = 4ir e(r')r' 2 dr' . We employ a fourth-order 
Runge-Kutta method together with a tabulated equation of state obtained from each of 
the holographic models. 

We first consider an EoS calculated in the D4/D8/D8 model. Physics of dense matter 
in Sakai-Sugimoto model [9] has been developed with/without the source term for baryon 
charge [14, 15]. One thing that might need to be improved with the Sakai-Sugimoto model 
for nuclear matter might be the absence of the scalar field which is responsible for the 
intermediate attraction of the nuclear force. In this work we adopt a simple EoS at zero 
temperature with a point-like baryon source given in [15]. Hereafter, we adopt the natural 
units. The regularized Helmholtz free energy is given by 

where rig is the baryon number density and 

K = 1 + Z 2 , a = 3.76 x 10 9 MeV 4 , b = 7.16 x 10~ 6 MeV -2 . 

Here they took A ~ 16.71 and M KK ~ 950 MeV [15]. We adopt the EoS in [15] to solve 
the TOV equation. In this case, however, we could not find a stable compact star, i.e., a 
star satisfying pressure-zero condition with a radius R, p(R) = 0, within a reasonable value 
of the radius. This means that the EoS from the D4/D8/D8 model may not support any 
stable compact stars or may support one whose radius is very large. This might be due 
to a deficit of attractive force in the D4/D8/D8 since the gravity alone may not balance 
against the pressure of the D4/D8/D8 EoS. In nuclear physics, the long range attractive 
force is governed by one-pion exchange, while the intermediate range ~ 1 f m attraction 
is dominated by a scalar field or two-pion exchange. In the D4/D8/D8 [9] a scalar that 
couples to baryon fields is missing. We may find a stable compact star if we try some other 
EoS based on Sakai-Sugimoto model other than the one given in [15]. However, the radius 
of the resulting star will be much larger than the one in nature and also the one obtained 
in the D4/D6 models given below. 

Now, we move on to the D4/D6 models. The EoS is given in Eq. (2.20) as 

p (p) = a(^n*l-[q]) (3.3) 

where p = b Q is the baryon number density and 

a ~ 4.45 x 10 8 MeV 4 , b ~ 1.07 x 10 6 MeV 3 . (3.4) 

Here we take A = 6 and M KK = 1.04 GeV in D4/D6 model [16, 17]. Since the value of 
Mkk is fixed uniquely by a meson mass, -Mkk ~ 1 GeV is a rather robust compared to 
the value of the 't Hooft coupling. With this choice of our model parameters, we plot the 
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(a) (b) 

Figure 1: (a): Relation between energy density and baryon number density, and (b): relation 
between pressure and baryon number density for the three holographic equations of state considered. 
Here we take A = 6 and M KK = 1.04 GeV. 



energy density and pressure with respect to the baryon number density in Fig. 1. For a 
given central baryon number density, a configuration is obtained by integrating Eq. (3.1) 
from the center, r = 0, to the point where the pressure vanishes, i.e., p(r = R) = 0, 
with R being the radius of the configuration. The mass of this configuration is given by 
m(r = R) = M. For low baryon number densities, i.e., at p < 3 x 10 -4 • po, the Baym- 
Pethick-Sutherland equation of state [18] is used. Here po is the normal nuclear density, and 
we use po = 1.3 x 10 6 MeV 3 . This equation of state is patched smoothly to the holographic 
equation of state at p ~ 3 x 10~ 4 • po MeV 3 . 

By solving the TOV equation, we obtain a mass-radius relation for each equation of 
state. In Fig. 2, we compare the mass-radius relations for three holographic equations 
of state with Nf = 1,2,3. As we increase the number of quark flavors we find that 
the maximum mass and corresponding radius becomes smaller. This is consistent with the 
simple picture that larger Nf favors a soft EoS, thereby smaller maximum mass and radius. 
Roughly speaking, the observed neutron stars have masses around 1.5M© and their radii 
are around 10 km. Compared to these observed neutron stars, what we obtained are more 
massive and larger. 

To see the A-dependence of our results, we take A = 17 and Mkk = 1-04 GeV. We 
show the energy density and pressure for each model with different values of A in Fig. 3. 
In Fig. 4, we collect mass-radius relations with A = 6, 17 and Mkk — 1-04 GeV. We find 
that with a bit larger value of the 't Hooft coupling, the mass and radius of our compact 
stars becomes smaller. For instance the maximum mass for Nf = 2 with A = 6 was about 
3.86Af Q and its radius was 45 km, which are now M = 1.36M Q and R = 17 km. 

So far we have focused on the mass-radius relation from holographic EoSs. Now we 
look into the holographic EoS in detail to see why our compact stars are different from 
observed neutron stars and to check if it satisfies constraints from thermodynamics or 
causality. In Fig. 5, we compare the pressure from the D4/D6 models with those from a 
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Figure 2: Comparison of mass-radius relations for three holographic equations of state with A = 6 
and M KK = 1.04 GeV. 




Figure 3: (a): Relation between energy density and baryon number density, and (b): relation 
between pressure and baryon number density for the three holographic equations of state considered. 



few conventional EoSs: APR [19], FPS [20], Sly4 [21]. Note that pressures from D4/D6 
models are quite larger than those from the conventional EoSs in most of the density region 
that are relevant to neutron stars. This might be the primary reason why our compact stars 
are quite large in radius compared to other conventional results. Large pressure implies 
large speed of sound. In order to check if the EoSs from D4/D6 models don't violate the 
causality, we calculate the square of the sound velocity, (? s = dp/de on the left panel of 
Fig. 6. Our result shows that the speed of sound is less than that for the fully relativistic 
Fermi gas limit, (? s = 1/3. We also calculated the adiabatic polytropic index 7 defined by 
7 = c?lnp/(dlne) on the right panel of Fig. 6. Note that the results with A = 17 shows 
lower speed of sound but higher adiabatic polytropic index compared to those with A = 6. 
This result can be understood from the relation 7 = c^(e/p) and the results in Fig. 3. 
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Figure 4: Comparison of mass-radius relations for three holographic equations of state. 
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Figure 5: Comparison of the pressure from D4/D6 models with those from a few typical EoSs. 



4. Summary 

We solved the TOV equation using the EoS calculated in holographic QCD. We first used 
the EoS from a dense D4/D8/D8 model. In this case, however, we could not find a stable 
compact star that satisfies pressure-zero condition p(R) = at a reasonable value of the 
radius R. This means that the EoS from the D4/D8/D8 model may not support any stable 
compact stars or may support one whose radius is very large. This might be due to the lack 
of attractive force mediated by a scalar field or two-pion exchange in the model. Next, we 
considered the D4/D6 models with different number of quark flavors, Nf = 1, 2, 3. Though 
the mass and radius of a holographic star are larger than those of normal neutron stars, 
we found stable compact stars in this case. To see why the compact stellar object obtained 
from the D4/D6 model is quite different from observed compact stars, we compared the 
pressure from our model with a few typical ones. We found that the pressure in our work 
is order of magnitude larger than that of conventional works. We also confirmed that the 
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Figure 6: (a) Sound velocity squared c 2 s , and (b) the adiabatic polytropic index 7, for the holo- 
graphic EoSs with respect to the normalized baryonic number density. 



EoSs considered in our work don't violate causality by considering the sound velocity square 
c 2 s and the adiabatic index 7. Since holographic QCD is based on the QCD symmetries 
which are believed to be restored at very high densities, extending the model to low density 
region may not be valid. In order to have a more realistic compact star, one may have to 
find a way of introducing proper attraction in holographic QCD or have to make a hybrid 
star in which low density EoS is replaced by the conventional one. 
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